Quantum phase transition in the sub-Ohmic spin-boson model 
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We study systematically the quantum phase transition (QPT) of the sub-Ohmic spin-boson model 
(SBM) in the full range of the power index s of the spectral density. Under the rotating-wave 
approximation (RWA), it is found that the usual delocalized-localized QPT is absent but a novel 
QPT occurs, which can be deduced analytically from whether a bound state between the spin and 
its reservoir exists or not. Surprisingly, this QPT has been missed in the literature for a long 
time. When the RWA is relaxed, the delocalized-localized QPT is recovered by using a unitary 
transformation plus perturbation method. Meanwhile, the novel QPT found under the RWA still 
exists but is located in the delocalized phase regime. In both cases the novel QPT causes a dynamical 
transition of the spin system from complete decoherence to decoherence suppression, which suggests 
a useful control way to the detrimental effects of the reservoir to the system. The result also implies 
that the coherent-incoherent transition studied extensively in the literature is caused essentially by 
the intrinsic QPT of the system. 
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I. INTRODUCTION 

Quantum phase transition (QPT), occurring at the ab- 
solute zero temperature, is induced by varying the exter- 
nal temperature-independent parameters, which results 
in a sudden qualitative change of the macroscopic prop- 
erties mapped from the eigen-spectrum of a quantum 
many-body system [1]. Recently, there are intense in- 
terests in QPT of a small quantum system coupled to a 
dissipative reservoir, e.g. the spin-boson model (SBM) 
and the Dicke model [2]. The SBM describes an effec- 
tive two-level system (TLS) interacting with a collection 
of harmonic oscillators acting as the bosonic reservoir. 
The effect of the reservoir is characterized by the spec- 
tral density J(w), which is phcnomcnologically described 
by the power law to the frequency of the reservoir, i.e. 
J(oj) oc oj s . The renewed interest in the SBM arises from 
the fact that the interaction of the TLS as a qubit with 
its reservoir leads to the decoherence which is a main ob- 
stacle to implement quantum information processing [3] . 
One desires that the exploration the QPT in the SBM 
can supply some insight to the decoherence control of the 
qubit system. A variety of approximate analytical and 
numerical methods, for example path integral method 
under non-interacting blip approximation [4] , variational 
method based on unitary transformation [5, 6], numerical 
renormalization group method [7] , quantum Monte Carlo 
method [8], and numerical diagonalization in coherent- 
state basis [9], have been developed. A consensus is that 
the SBM shows a QPT from derealization to localization 
with the increase of the dimensionless coupling constant 
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in the Ohmic case (s — 1), as a consequence of the com- 
petition between the internal tunneling effect of the TLS 
and the external dissipation effect of the reservoir. 

The QPT of the SBM in sub-Ohmic case (0 < s < 1) 
has attracted much attention recently. One of the mo- 
tivations originates from the idea that the sub-Ohmic 
reservoir can be used to model l/f noise [10], which is 
the main source of decoherence in solid state systems, 
e.g. quantum dots [11, 12] and superconductor qubit sys- 
tems [13, 14]. Another motivation is that the sub-Ohmic 
SBM is more involved than the Ohmic case. Different 
methods to the sub-Ohmic SBM can even not lead to a 
qualitatively consistent result. The path integral method 
under non-interacting blip approximation predicts that 
the QPT from dclocalization to localization is absent for 
the sub-Ohmic SBM [4]. The numerical renormalization 
group method confirms the occurrence of the QPT in 
the full range < s < 1, while the breakdown of the 
quantum-to-classical mapping for < s < 1/2 [7], which 
means the failure of the classical mean-field description 
to the QPT. However, the quantum Monte Carlo method 
and the numerical diagonalization in coherent-state ba- 
sis predict the presence of the QPT, the well-defined 
quantum-to-classical mapping, and the classical critical 
exponents to the sub-Ohmic SBM [8, 9]. 

While the physics in the localized regime of the SBM is 
trivial, the delocalized regime exhibits many interesting 
phenomena. For example, it was shown that the dynam- 
ics of the spin in the delocalized regime changes from 
the damped coherent oscillation to incoherent relaxation 
with the increase of the coupling strength for both of the 
Ohmic [4] and the sub-Ohmic [15-17] SBM. This dynam- 
ical phenomenon is named as the coherent-incoherent 
transition [4]. What is the physical reason of this dy- 
namical transition? As a dynamical behavior, it does 
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not suffice to deduce the occurrence of the intrinsic QPT 
of the model because it is initial-state-dependent. By 
evaluating the ground-state energy of the Ohmic SBM, 
it was shown in Ref. [18] that the coherent-incoherent 
transition is actually the dynamical consequence of an- 
other intrinsic QPT occurred in the dclocalized regime. 

In the present work, we study the ground state energy 
and the dynamics of the TLS of the sub-Ohmic SBM 
both with and without the rotating-wave approximation 
(RWA). Under the RWA, we find analytically that a novel 
QPT occurs in the full range of the sub-Ohmic spectrum, 
while the delocalized-localized QPT is absent. This QPT 
separates the dynamics of the TLS into complete deco- 
herence and decoherence suppression. When the RWA 
is relaxed, using the perturbation approach based on a 
unitary transformation, which has been successfully used 
to capture the delocalized-localized QPT for the Ohmic 
[5, 6] and the sub-Ohmic [17] SBM, we find that the 
novel QPT happened explicitly in the conventional dclo- 
calized phase regime still exists. It implies that besides 
the conventional delocalized-localized QPT, there is an- 
other QPT in the SBM. The compatibility of the novel 
QPT with the coherent-incoherent transition makes us 
conjecture that the coherent-incoherent transition occur- 
ring in the delocalized regime is actually caused by an in- 
trinsic QPT. Our analytical formulation provides a clear 
physical picture of this QPT and a unified description of 
the QPT in the sub-Ohmic SBM. 

This paper is organized as follows. In Sec. II, the SBM 
and its simplification under the RWA arc introduced. To 
illustrate the new kind of QPT, we explore in Sec. Ill 
the QPT in the SBM under the RWA by examining the 
formation of a bound state. In Sec. IV, the QPT in 
the sub-Ohmic SBM is investigated by means of the per- 
turbation approach based on a unitary transformation 
[5, 6]. A novel QPT is found in the usual delocalized 
phase regime. Sec. V is devoted to the dynamical con- 
sequence of such QPT. Finally, a brief discussion and 
summary are given in Sec. VI. 



II. THE MODEL 

The SBM, which relates to a variety of physical and 
chemical processes [19], describes the tunneling of a quan- 
tum particle in a double well potential under the influence 
of a bosonic reservoir. Its Hamiltonian reads 

H = - j&* + £ uj k b{b k + £ f& z (b k + V k ), (1) 

k k 

where e is the difference between the two energy lev- 
els in the wells, A is the tunneling amplitude between 
the two wells, b\ and b k are, respectively, the creation 
and annihilation operators of k-th mode of the reser- 
voir with frequency uj k . The coupling strength between 
the particle and its reservoir is denoted by g k , which 
is further characterized by the spectral density J(uj) = 



7r E/clfl'fc| S(co — u} k ). In the continuum limit the spectral 
density may have the form 

J(cj) = 2Trauj 1 c - s uj s e(Lu c - cj), (2) 

where a is a dimensionless coupling constant, uj c is a 
cutoff frequency, and Q(x) is the usual step function. 
The reservoir is classified as Ohmic when s = 1, sub- 
Ohmic when < s < 1, and super-Ohmic when s > 1 
[4]. In spite of the simplicity of its formulation, the SBM 
does not admit an exact solution in a closed analytical 
form and one often resorts to numerical simulations or 
various approximations for its analysis. Under a unitary 
transformation Ui = cxp(— iiray/A), one can prove that 
Eq. (1) is equivalent to 

H z = + + J2u k blb k + ^2^-a x (b k + b{), (3) 

k k 

which corresponds to a it/ 4- rotation around a y to SBM 
(1). In the following, we assume e = 0, which corresponds 
to the case of the un-biased double potential wells where 
there has no energy difference between the two wells. 

The interaction in Eq. (3) contains the counter rotat- 
ing terms, b\,<j-\. and b k a~. A widely used approximation 
in quantum optics and quantum information communi- 
ties is the RWA, which is applicable in the weak coupling 
limit. Then Eq. (3) is reduced to 

k k 

Eq. (4) is analytically solvable because the total excita- 
tion number N — ^Z k b\b k + &+<J- of the whole system 
is conserved. 

III. QPT UNDER THE RWA 

The Hamiltonian (4) under the RWA is widely used 
to characterize decoherence of a qubit in quantum optics 
and quantum information communities. It is found that 
the decoherence dynamics of the TLS in this situation 
exhibits rich behaviors in different spectral densities of 
the reservoir, such as the exponential decay under the 
Born-Markovian approximation [20], the oscillatory de- 
cay in a lossy cavity with Lorentzian spectral density [21], 
even the decoherence suppression in engineered reservoir 
with photonic band gap structure [22-24] . Further stud- 
ies show that whether a bound state of the global sys- 
tem is formed or not essentially plays a key role in such 
rich dynamical behaviors [25-27]. Experimentally, the 
bound-state induced decoherence suppression has been 
observed [28-30]. In this section, starting from Eq. (4), 
we investigate analytically and numerically QPT via ex- 
amining the formation of a bound state between the TLS 
and its reservoir. We reveal that the formation of the 
bound state actually implies a QPT of the system. 
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A. Analytical results 

Since N is conserved, the Hilbert space is split into 
the direct sum of the subspaccs with definite quantum 
number N. In this situation one can naively deem that 
the eigcnstate \ipo) = |— ,{0fe}}, a tensor product of the 
respective ground states of the two subsystems in zero- 
excitation subspace with eigenvalue Eq = —A/2, is the 
ground state of the whole system. Is this always true? 
To verify this, let us examine the eigen solution of //tota 
in the single-excitation subspace, which can be expanded 
as \<px) = c o \ + ,{0 k }) + Y / ^Lo c kh> 1 k)- From the eigen- 
cquation governed by Eq. (4) we can obtain a transcen- 
dental equation of E\ 



J(U>) 



-dui = Ei. 



(5) 



A bound state is an eigenstate with real (discrete) eigen- 
value in a quantum many-body system. So if Eq. (5) 
has real root, we can claim that the system possesses 
a bound state [26, 31]. We can easily find that y(E\) 
decreases monotonically with the increase of E\ in the 
regime of E\ < — ~. Therefore if the condition 



tf(-2-)<-2. 



(6) 



is satisfied, y(E\) always has one and only one intersec- 
tion with the function on the right-hand side of Eq. (5). 
This root just corresponds to the eigenvalue of the formed 
bound state in the Hilbert space of the system plus its 
reservoir. On the other hand, in the regime of E\ > — A , 
we can see that y(E{) is divergent, which means that 
no real root Ei can make Eq. (5) well-defined. Conse- 
quently, Eq. (5) does not have real root to support the 
existence of a further bound state in this regime. It is 
noted that Eq. (5) may possess complex root. Physi- 
cally, this means that the corresponding eigenstate ex- 
periences decay contributed from the imaginary part of 
the eigenvalue during the time evolution, which causes 
the excited-state population approaching zero asymptot- 
ically and the decoherence of the TLS. However, to the 
bound state, the excited-state population is constant in 
time. This means that the formation of bound state can 
result in decoherence suppression. 

With the formation of the bound state, the ground 
state is changed from \ipo) to the bound state \<pi), be- 
cause Ei < Eq. Wc can verify that the eigenstates of 
-Hrwa in the subspaces N >2 actually have larger eigen- 
values than Ei, which is shown in Appendix A. This im- 
plies that the higher-boson states may not become the 
ground state. One may also observe that the two states 
arc orthogonal, i.e. (—,{0k}\<fi) = 0. Therefore, the 
energy-level crossing accompanying with the formation 
of bound state signals clearly that the system undergoes 
a QPT. From the criterion (6), it is readily to evaluate 
that the QPT happens at the critical point 

_2sA 




FIG. 1. (Color online) Ground-state energy E g (a) and its first 
derivative E' g — (b) as a function of the coupling constant 
a and power index s of the spectrum. The parameter used 
here is A = 0.1oj c . According to Eq. (7), the first-order QPT 
occurs at a c RWA = 0.2s, which has been confirmed by the 
discontinuousness of the first derivative of the ground-state 
energy. 
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FIG. 2. (Color online) Ground state fidelity (a) and entan- 
glement entropy (b) as a function of the coupling constant a 
and power index s of the spectrum. Here 8a — 0.0005, other 
parameters used here are the same as Fig. 1. The singularity 
near the critical point a c RWA of forming bound state shows 
existence of QPT in this model. 



for the spectral density (2). 



B. Numerical results 

To verify the existence of QPT in the model quantita- 
tively, in the following we study numerically the ground- 
state energy and its derivative, fidelity and entanglement 
entropy of the ground state near the critical point with 
the change of the coupling constant of the system. 

At zero temperature, the non-analyticity of the 
ground-state energy is directly connected to the QPT. 
The first (or n-th) order QPT is characterized by the dis- 
continuity in the first (or n-th) derivative of the ground 
state energy. In Fig. 1, we plot the ground state en- 
ergy and its first derivative. It can be seen that the first 
derivative is discontinuous at the critical point (7), which 
means that it is a first-order QPT. 

The QPT can be further verified by the fidelity F and 
entanglement entropy £ between the TLS and the reser- 
voir of the ground state. The ground-state fidelity is 
defined as the overlap of two ground states correspond- 
ing to two slightly different control parameters [32], i.e. 
F = (tpg(a)\tp g (tp + 8a)). The entanglement entropy can 
be obtain by calculating the entropy of the reduced den- 
sity matrix of the TLS after tracing out the reservoir de- 
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grees of freedom. In Fig. 2(a), we plot F near the critical 
point. The singularity in the plot evidences clearly the 
existence of QPT in this model. Because of the totally 
orthogonal property of ground state, the fidelity com- 
pletely drops to zero at the critical point. In Fig. 2(b), 
we plot £ of the ground state. Near the critical point, 
wc find a sudden birth of the ground-state entanglement, 
which can be seen as a result of changing of ground-state 
structure. This discontinuousness in the ground-state en- 
tanglement entropy also evidences the existence of QPT. 



IV. QPT WITHOUT THE RWA 

In the following, using the perturbation approach 
based on unitary transformation [5, 6], we study whether 
or not the QPT still occurs when the RWA is relaxed. 



A. Analytical results 



A unitary transformation J7 2 = exp[J^ fc — 



bk)&x] can recast Eq. (3) into 



k k 

A 



E5fc(l - 6c) ,t , tu» 
^ {bk + K)ir x 



[a z cosh x — iffy sinh x] + C, 



(8) 



where C = £ fc - 2), * = E k - h), and 

6c are to be determined. Wc can further separate H' into 

H' = H' Q + V', (9) 

where 

■qA 



Hl ) = ^-a z +Y,^kb{b k + C, 



A 



(b k + b\)fr x - i—Oy sinh x 



A 



+ y^( cosh X - rj), 



(10) 



with 



r, = ({0 fe }|coshx|{0 fc }> =cx P [-^0]. (11 



We can see from Eq. (9) that, to the zero-order ap- 
proximation, the spin-boson interactions can be elimi- 
nated to generate an effective nonintcracting Hamilto- 
nian characterized by a renormalized tunneling ampli- 
tude A e ff = ijA with rj as renormalized factor. 

With the separation of Eq. (9), we can readily calcu- 
late the Bogoliubov-Peierls bound on the free energy Fb 
of the system [33]. The free energy F of the system is 
related to Fb by F < Fb with 



Fr, 



-p-'teTrexpi-pHti + iV')^, (12) 



where (3 = j^p, (')it' denotes the thermal expectation 
value calculated with respect to H[ ): and the trace is cal- 
culated using the eigenstates of H' . It is easy to find 
(V')h = 0. Therefore, 



ln[2cosh£|2] 

is 



c. 



(13) 



The parameters are determined by minimizing Fb 
with respect to that is = 0. Wc find in our 

zero-temperature case (i.e. (3 — > 00) 



Lo k + r]A 



(14) 



By now, the parameters as well as the renormalized 
factor rj have been determined. The renormalized factor 77 
has been used successfully to characterize the delocalized- 
localized QPT in the SBM [5, 6, 17]. If the tunneling 
amplitude is renormalized to zero, the system is in the 
localized phase and the dynamics is trivial. In contrast, if 
the renomalizcd tunneling amplitude is nonzero, then the 
system is in the delocalized phase, which displays some 
interesting dynamical behaviors such as damped coherent 
oscillation and incoherent relaxation. 

Focusing on the delocalized phase regime, we further 
separate the first-order perturbation term from V' = 



H[ + H' 2 with 



z {b k + b k )a x - i—(J y X 



^v k (b k a + + &j>_), 



H'2 = y<Mcoshx- rj) - i—a y {smhx - rjx), (15) 

where v k = r]Ag k ^ k /ui k and a± = {a x ± ia y )/2. Com- 
bined with Eq. (9), we arrive at the transformed SBM 
as H' = Hq + H [ + H' 2 , where H' collects all the renor- 
malized non-interacting terms, H[ collects all the first- 
order perturbation terms, and H 2 collects all the higher- 
order perturbation ones. It means that by appropriately 
choosing the unitary transformation, we can separate au- 
tomatically the excitation-conservative transition terms 
Hq and H[ from the multi-boson-excitation and nondi- 
agonal transition H 2 . It has been proved that in zero- 
temperature and weak coupling (i.e. the delocalized) 
regimes the higher-order perturbation terms H 2 can be 
neglected [6] . Then the transformed Hamiltonian has the 

-ffeff 



form W &H' Q + H[ 



H cfi = ^-a z + J2[ukb{h + v k {b\a- + b k a+)] + ($6) 
k 

which shares the formal similarity with the rotating-wave 
approximate Hamiltonian (4). 

With the similar procedure as in Sec. Ill and neglect- 
ing temporarily the constant term C in Eq. (16), we 
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can determine that a bound state \(p\) = d a |+,{0}) 
Tlik^k \— j {l}fe) with the eigenvalue Ei satisfying 



y(Ei) 



77A 

~2~ 



E 



= £1 



(17) 



can be formed for H c g . This equation permits a real root 
in the regime E\ < —77 A/2 if and only if y{— T]A/2) < 
—j]A/2. Accompanying with the formation of a bound 
state, the ground state is changed from \ip' Q ) = |— , 0^) to 
Recovering back the neglected term C, we get the 
ground-state energy as 



C, a < a r 

2 

Ei — C, a > a c 



(18) 



where the critical point a c can be determined by solving 
equation y(— 77A/2) = — 77A/2. Physically, such a sudden 
change of the ground-state structure signals clearly the 
occurrence of QPT in the system. 

It is noted that the neglected higher-order perturbation 
term H' 2 gives no contribution to the QPT because it is 
zero, i.e. ((p'^H^tp'j) = = 0,1), in the two eigen 

bases. It means that the neglected term H' 2 has no impact 
on such level-crossing-causcd QPT. This in turn validates 
our approximation. 



B. Numerical results 




FIG. 3. (Color online) The conventional delocalized-localized 
QPT characterized by the renormalized factor rj as a function 
of the coupling constant a and the power index s of the spec- 
trum. The red solid line depicts the critical point. A = 0.1ui c 
has been used in the numerical calculation. 

For completeness, we firstly recover the conventional 
delocalized-localized QPT in the full range of the sub- 
Ohmic spectral density. In Fig. 3, we plot the numeri- 
cal results on this conventional QPT characterized by 77, 
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FIG. 4. (Color online) (a): Ground-state energy E g and (b): 
its first derivative E' g = as a function of a and A, and (c): 
the phase diagram when s = 0.8. The red dashed line and 
the blue solid line mark the critical points of the conventional 
delocalized-localized QPT and the novel bound state QPT, 
respectively. 



which can be calculated by solving Eqs. (11) and (14) 
self-consistently. We can see that the system is in the 
delocalizcd phase regime when a is small, where 77 takes 
a finite value. With the increase of a, 77 drops suddenly 
to zero and the system enters into the localized phase 
regime. Such delocalized-localized QPT is present in the 
whole range of the power index s of the sub-Ohmic spec- 
tral density. This is coincident with the results under the 
quantum Monte Carlo method and the numerical diago- 
nalization method [8, 9]. The critical point of this QPT 
tends to a c — 1 with the increase of s to Ohmic case. 
This is consistent with the well-known results that the 
delocalized-localized phase transition occurs at a c = 1 
in the small A limit for the Ohmic SBM [4] . 

Focusing explicitly on the conventional delocalizcd 
phase regime, where r\ takes a finite value, we now study 
the ground state energy of the sub-Ohmic SBM. From the 
analysis above, we know that the formation of a bound 
state in the single-excitation subspace of H c fi causes the 
level crossing of the ground state, which triggers a novel 
QPT to the model. To verify this, taking s = 0.8 as 
an example, we plot in Fig. 4(a) and (b), respectively, 
the ground state energy and its first derivative to the 
coupling constant according to Eq. (18). We can see 
that Eg is continuous, but -j- 3 - shows a discontinuity at 
the critical point a c where the bound state is formed. 
It manifests clearly that there is another QPT existing 
in the delocalized phase regime. In Fig. 4(c) the phase 
diagram is depicted. It indicates that besides the conven- 
tional delocalized-localized QPT with the critical point 
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FIG. 6. (Color online) The nonequilibrium dynamics of P z (t) 
of the SBM with the RWA under the initial condition Eq. 
(19). Here the parameters s = 0.8 and A = 0.1o; c have been 
used. Consequently, the critical point is ac.RWA = 0.16. 



FIG. 5. (Color online) (a): Ground-state energy E g and (b): 
its first derivative E' g = as a function of a and A, and (c): 
the phase diagram when s = 0.4. The red dashed line and 
the blue solid line mark the critical points of the conventional 
delocalized-localized QPT and the novel bound state QPT, 
respectively. 

marked by the red dashed line, there is another novel 
QPT with the critical point marked by the blue solid 
line occurring in the delocalizcd phase regime. To check 
if such novel QPT is also present in < s < 1/2, where 
the quantum-classical mapping was proven to be break- 
down by the numerical renormalization group method 
[7], we take s = 0.4 as an example and plot explicitly 
the ground state energy, its derivative, and the phase di- 
agram in Fig. 5. We can see although the delocalized 
phase regime much shrinks with the decrease of s, it still 
separates into two further regimes with and without the 
bound state. It shows that the novel QPT also shows up 
when < s < 1/2. 



V. DYNAMICAL CONSEQUENCES OF THE 
QPT 

To evaluate the consequences of the novel QPT on the 
non-equilibrium dynamics of the spin system, we study 
next the evolution of the spin system under certain ini- 
tial condition. It has been found that the spin dynamics 
shows the so-called coherent-incoherent transition [4] in 
the delocalized phase regime when the initial state of the 
reservoir is in vacuum [15, 16]. Some theoretical [22-27] 
and experimental [28-30] works also show that the dy- 
namics of the SBM under the RWA exhibits rich dynam- 
ics, from complete decoherence to dccohcrcnce suppres- 
sion. We argue that both of the two different dynamical 
behaviors are essentially caused by a same reason, i.e. 



the novel QPT wc revealed in Sec. Ill with the RWA 
and the in Sec. IV without the RWA. To verify this, we 
study in the following the dynamics of the SBM both 
with and without the RWA. 

Firstly, wc study the nonequilibrium dynamics of the 
spin under the RWA. Assuming initially the whole system 
is in 

l*i(0)> = |+,{0 fc }>, (19) 

which, in the H z , i.e. Eq. (3), representation, takes 
the form as |*,(0)) = #i|*i(0)) = |+.„{0 fc }) with |+ x ) 
satisfying a x \+ x ) = \+ x ). The time evolution of |5 , z (0)) 
is governed by Eq. (4) under the RWA. Therefore, the 
time-dependent solution can be expanded as 

\V z (t))=e^{-L\-,0 k )+<t)\+,0k)+Y;Mt)\-,l k )}- 

(20) 

From the Schrodingcr equation, we can get the probabil- 
ity amplitude c(t) satisfying 

c{t)+iAc(t)+ f f(t - t)c{t)cLt = 0, (21) 
Jo 

where the kernel function f(t — r) = 
Jo °° J( UJ ) e ~ tlJ< ' t ~ T ' > duj and the initial condition 
c(0) = l/v2- With this result at hands, we can calculate 
P z = (y(0)\e iAt & z e- i6t \i()(0)) under the RWA as 

P z it) = (tt,(t)|*x|tf,(t)} = \/2Re[c(i)]. (22) 

In Fig. 6 we plot P z {t) of Eq. (22) in different cou- 
pling constants. According to Eq. (7), the QPT occurs 
at a c RWA = 0.16. We can see that the dynamics reduces 
the oscillatory damping to zero when a < 0.16. It is un- 
derstandable from the fact that the bound state in this 
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region is absent and all the quantum coherence decays to 
zero. We call the character of the dynamics in this region 
as the complete decoherence. When a = 0.16, it is very 
interesting to find that the coherence does not decay to 
zero and a finite quantum coherence is preserved in the 
steady state. In this situation, the bound state \ipi) with 
the eigenvalue being just E\ = —A/2, which is equal to 
Eq, is formed. As a stationary state, the quantum coher- 
ence contributed from the component to Eq. (19) 
does not change during the time evolution. Therefore, 
we get a finite asymptotical P z (t). With the further in- 
crease of a, the dynamics shows lossless oscillation. The 
bound state with smaller eigenvalue than Eq is present. 
In this region, the components of |cpo) and \(fi) in Eq. 
(19) have different time dependence. The difference of 
the two eigenvalues, i.e. Eq — E±, contributes to the fre- 
quency of this lossless oscillation. A larger a induces a 
smaller E\ and a more faster oscillation of P z (t) . We call 
the character of the dynamics in the region a > a c RWA 
as the decoherence suppression. 

Next, we study the dynamics when the RWA is relaxed. 
To make the dynamics manifest the effect of the formed 
bound state exclusively, we choose the initial state as 

|* 2 (0)) = |+> ® tflKO*}). (23) 

It is noted that this state is different from Eq. (19), un- 
der which it has been shown that the dynamics of the 
SBM without the RWA exhibits the coherent-incoherent 
transition. The merit of choosing this state as the initial 
state is that it takes the form as |^(0)) = (7 2 f/i|* 2 (0)) = 
|+x,{0fe}) in the H e s representation. Therefore, besides 
the zero-excitation, only the single-excitation subspace 
where the bound state is formed is involved in the dy- 
namics. In the same manner as the above RWA case, we 
can calculate P z (t) = \/2Re[h(t)], where h(t) satisfying 

h(t) + ir)Ah(t) + [ f'(t - r)/i(r)dr = 0, (24) 
Jo 

with the initial condition being h(0) = and the kernel 

function fit — r) = j- / °° J'(cj)e~ tL ^ t_r )cLj connecting 
to the renormalized spectral density J'(w) = ^ fc v^5{ui — 

Fig. 7 portrays P z {t) under the initial condition (23) 
when the RWA is relaxed. We can see that the simi- 
lar behavior as Fig. 6 is present. When a < a c , the 
bound state is absent. Therefore, the dynamics shows 
complete decoherence with the quantum coherence de- 
caying to zero. When a = a c , a finite steady P z (t) can 
be obtained asymptotically due to the double degeneracy 
of |v?o) an( i the formed bound state \<p[). When a > a c , 
P z (t) shows decoherence suppression with the quantum 
coherence tending to lossless oscillation, where the fre- 
quency of this oscillation is determined by the difference 
between the two eigenvalues of \ip' Q ) and \<p[)- 

From the analysis above, we can see that the novel 
QPT induced by the formation of the bound state has 



t 1 1 1 r 




10 20 30 40 

cot/100 

c 

FIG. 7. (Color online) The nonequilibrium dynamics of P z (t) 
of the SBM without the RWA under the initial condition Eq. 
(23). Here the parameters s = 0.8 and A = 0.1ai c have been 
used. The critical point can be evaluated numerically at ac = 
0.258. 

a profound consequence to the nonequilibrium dynamics 
of the TLS in the conventional dclocalized phase regime. 
It induces a dynamical transition from complete deco- 
herence to decoherence suppression for the initial states 
where only the single-excitation subspace is involved. On 
the other hand, a widely studied case is the nonequi- 
librium dynamics of the TLS when the reservoir is ini- 
tially in a vacuum state, where the spin dynamics shows 
a transition from damped coherent oscillation to inco- 
herent relaxation, i.e. the so-called coherent-incoherent 
transition, in the delocalized-phase-regime. Therefore, it 
is reasonable to conjecture that the coherent-incoherent 
transition is also a dynamical consequence of this novel 
QPT on the initial vacuum state of the reservoir. This 
has been proved analytically for the Ohmic spectral den- 
sity in Ref. [18] that a c of the QPT matches well with 
the point of the coherent-incoherent transition. Thus, 
we can conclude that both of the transitions from com- 
plete decoherence to decoherence suppression and from 
damped coherent oscillation to incoherent relaxation are 
actually the different dynamical consequences of the same 
intrinsic QPT on different initial states. 



VI. CONCLUSION 

In conclusion, we have studied the QPT of sub-Ohmic 
SBM both with and without RWA. When the RWA is 
used, we reveal that a novel QPT induced by the forma- 
tion of a bound state in the single-excitation subspace 
occurs, while the delocalized-localized QPT is absent. 
When the RWA is relaxed, we first obtain that the con- 
ventional delocalized-localized QPT occurs in the whole 
range of the sub-Ohmic spectrum, which is consistent 
with the results under the quantum Monte Carlo method 
and the numerical diagonalization method [8, 9]. Then 



8 



using the perturbation approach to neglect the high-order 
interaction terms in unitarily transformed Hamiltonian, 
we show that the bound-statc-induccd QPT still exists in 
delocalized phase regime. The approximation is justified 
by that we are working in the weak coupling (i.e. the de- 
localized phase) regime and in the zero-temperature case, 
where the high-order excitations are negligible. On the 
other hand, one also can verify that the neglected terms 
give no contribution to the QPT we obtained, which in 
turn validates our approximation. We also have stud- 
ied the dynamical consequences of the novel QPT. It is 
shown that the QPT causes a dynamical transition from 
complete decoherence to decoherence suppression to the 
initial state in which only the single-excitation subspacc 
is involved. This result is compatible to the coherent- 
incoherent transition which happens to the state where 
the reservoir is initially in vacuum. It conjectures that 
the coherent-incoherent transition reported in the liter- 
ature is essentially caused by the intrinsic QPT of the 
SBM. 

As a final remark, the results obtained in the present 
work is only in the un-biascd SBM, i.e. e = 0. It is 
also interesting to generalize our discussions to the biased 
case. This work is in progress. 
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Appendix A: The bound state with higher-excitation 

The cigenstatc of H c s with N = 2 can be expanded as 
1*2) = J2(d k \+,l k )+e k \-,2 k } + fkM- 

k k'^k 

(Al) 



From the eigen equation H c ff\(f 2 ) = ^2 1 V2 ) , we obtain 
the following coupled equations 



End 



2«fc 



77A 



u k )d k + V2v k e k + ^2 v k>fk,k>, (A2) 



E 2 e k = V2v k d k + (2u k - ^)e fc , (A3) 
77A 

E 2 fk,k' = v k >d k + (uj k +uj k > —)/*,#• (A4) 

Substituting the solutions of Eqs. (A3) and (A4) in term 
of d k into Eq. (A2), we have 



77A 

E 2 = — + u} k 



E 2 -(-^ + 2oj k ) 



k' 



E 2 -{-^+uj k +u Jk ,y 



(A5) 



Generally, as the density of states per unit frequency is 
proportional to M, where M is the total number of bath 
oscillators, the microscopic coupling strengths have to 
scale as M" 1 / 2 in order to ensure that the spectral den- 
sity is well defined in the thermodynamic limit. There- 
fore, in the thermodynamic limit, M — > 00, v k itself is 
small [34], then Eq. (A5) reduces to 



E2 = -y- + w fc + 



k i E 2 - (-^ + oj k + Wfc/) 



Setting E 2 — E' 2 + uj k , we have 



^u k ,-{E> 2 + rfy 



(A6) 



(A7) 



which is the same as Eq. (6) in the main text. It means 

E' 2 =E l <E 2 . (A8) 

Therefore, even a bound state in the higher-excitation 
subspace can be formed, its eigenvalue is still larger than 
E\. This means that the higher-excitation bound state 
is not possible to be the ground state. 
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